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Central spin decoherence is useful for detecting many-body physics in environments and moreover, the spin 
echo control can remove the effects of static thermal fluctuations so that the quantum fluctuations are 
revealed. The central spin decoherence approach, however, is feasible only in some special configurations 
and often requires uniform coupling between the central spin and individual spins in the baths, which are 
very challenging in experiments. Here, by making analogue between central spin decoherence and 
depolarization of photons, we propose a scheme of Faraday rotation echo spectroscopy (FRES) for studying 
quantum fluctuations in interacting spin systems. The echo control of the photon polarization is realized by 
flipping the polarization with a birefringence crystal. The FRES, similar to spin echo in magnetic resonance 
spectroscopy, can suppress the effects of the static magnetic fluctuations and therefore reveal dynamical 
magnetic fluctuations. We apply the scheme to a rare-earth compound LiHoF4 and calculate the echo signal, 
which is related to the quantum fluctuations of the system. We observe enhanced signals at the phase 
boundary. The FRES should be useful for studying quantum fluctuations in a broad range of spin systems, 
including cold atoms, quantum dots, solid-state impurities, and transparent magnetic materials. 

When a central spin is coupled to a spin bath, the quantum coherence of the central spin would be lost 
due to the noise from the bath^'^'^. On the one hand, the central spin decoherence is an important issue 
in quantum coherence based technologies such as quantum computing^'^ and magnet ometry^"^. On 
the other hand, since the noise is caused either by thermal fluctuations or by elementary excitations in the bath, 
the central spin decoherence is a useful probe of many-body physics in the bath. A number of interesting effects 
have been discovered by viewing the central spin decoherence as a probe of many-body physics^°"^^. 

Coupling a single spin to a spin bath, however, is not a trivial task in experiments and is feasible only in certain 
specially designed systems. Moreover, such coupling between the central spin and individual bath spins is often 
required to be uniform^°"^^, which puts additional constraints on experiments. A potential solution is motivated 
by the spin noise spectroscopy^^"^^, where the polarization of photons plays the role of a central spin and the 
depolarization of the photons due to the spin noise resembles the central spin decoherence. The coupling between 
photons and spins has been well established for many systems including atoms^^, quantum dots^^'^^, and trans- 
parent magnetic materials^". A famous example is the Faraday rotation Due to slow spatial variation of laser 
pulses, the coupling can be easily made almost uniform for all spins interacting with the photons. Therefore, in 
this paper, we propose a scheme based on Faraday rotation and photon depolarization to study many-body 
physics of spins systems, in lieu of the central spin decoherence method^°"^^. 

A particularly useful technique in central spin decoherence method, which is not yet available in the optical 
spin noise spectroscopy, is the spin echo^^, or more generally, dynamical decoupling controP^' In spin echo or 
dynamical decoupling, the central spin is flipped and effectively senses the noises from opposite directions before 
and after the flip control. Therefore, if the noise is static (or slow enough), the effect would be cancelled when the 
time the central spin evolves after the flip equals that before the flip and the central spin coherence will recover at 
that particular time, resembling an echo. Since the static noises usually result from the classical configurations of 
the laboratories and from the thermal fluctuations in the baths, the central spin decoherence under the spin echo 
control is particularly sensitive to the dynamical fluctuations in the bath, which have quantum origins. Therefore, 
it has been shown before that even at a high temperature, the central spin decoherence, when it is under echo 
control, can still single out the quantum fluctuations^^. To incorporate the echo control in the optical noise 
spectroscopy, we design an optical set up in which the photon polarization can be flipped by a birefringence 
crystal and therefore the analogue to spin echo can be realized in the photon depolarization measurement. We test 
our scheme on a rare-earth compound LiHoF4 which is an experimental realization of Ising magnets^^'^^ and has 
sizable Faraday rotation effect on optical probes^°. 
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Results 

Fluctuations of Faraday rotation angle. A central spin- 1/2 (|t) + 
1 1)) under an external field b{t) along the z- axis accumulates a phase 

shift such as (e'^|t) + with 6= b(t)dt. If the external field is 

Jo 

random (due to interaction with fluctuating spins in the bath), the 
coherence of the spin will be lost. 

In analogue to this spin decoherence, we consider depolarization 
of a linearly polarized laser pulse. Being a superposition of two cir- 
cular polarizations (cr^ + cr"), the laser pulse can play the role of a 
central spin. When the linearly polarized laser pulse propagates 
through a magnetic sample that contains fluctuating spins (J/) [see 
Figure 1 (a)], the two circular components get different phase shifts 
{e'^cT^ + e~'^a~), which produce a rotation angle 0 of the linear 
polarization, known as the Faraday rotation (FR). The FR angle is 
proportional to the magnetization along the direction of propagation 
and the thickness (L5) of the sample. Here we consider the case that 
the pulse length is much longer than the thickness of the transparent 
sample. The laser, therefore, can be considered to be interacting with 
all the spins within the cross sectional area (A) simultaneously. The 
FR angle gained by a linearly polarized laser after it travels through 
the sample along the z direction is 



0(tp) = at-'Ls 



Mz(t)dt, 



(1) 



where tp is the duration of the laser pulse, Mz(t) = N~^ ^ Jf(t) is the 

/-I 

magnetization, N = pALs is the number of spins located within the 



(a) 



(b) 




fhlX 



interaction region (p being the density of spins), and a is a coupling 
coefficient, which depends on the laser frequency and the material 
parameters (see Methods). 

The fluctuation of the magnetization causes a random phase shift 
of the circular polarized components of the laser, leading to the 
depolarization of the laser, which is similar to the decoherence of a 
spin- 1/2 under a random field. The fluctuation of the FR angle is 
determined by the fluctuation of the magnetization Czih — t^ = 
{AMz(ti)AMz(t2)) as (see Methods) 



dtidtiCziti — ^2)? 



(2) 



where A^fr = ^fr (^fr) and AM^ = - (M^). The symbol (...) 
denotes the ensemble average, that is, {M)= y^P„(i/^^|M|i/^^), 

n 

where and are, respectively, the eigenstate and eigenenergy 
of the spin system, and Pn = e~^^" / e~^^'^ is the probability dis- 

m 

tribution, with P being the inverse temperature. The Fourier trans- 
form of the correlation function Q(0 gives the noise spectrum of the 

spin system S(co) = | Q(0 exp (iwt)dt. 

The noise spectrum can be understood as caused by two mechan- 
isms. One part of the fluctuations is caused by the transitions between 
different energy eigenstates. The corresponding spectrum is 

Sq((j)) = 2k^^^ 3[(o — (Eri—Ern)]Pn\{^n\^\^rn)\^^ which contains 

non-zero frequency components. Therefore, it is dynamical and 
quantum. Another part of the fluctuations is originated from the prob- 
ability distribution P„ at finite temperature and has only the zero 

frequency component: St\,(co)=2nS(w)^2Pn{^n\^\^n)^ ~ (^)^^ ' 

which is the static thermal fluctuation and vanishes at zero temper- 
ature. The effects of the static thermal fluctuations can be removed by 
spin echo^^, or by the Faraday rotation echo as studied in this paper. 
Therefore the features of the quantum fluctuations would be revealed 
by the echo methods. 



(c) 



(d) 



Birefringence crystal 




f/ec, 



Figure 1 | Schematic scheme of Faraday rotation echo, (a) A linearly 
polarized laser pulse passes through a sample containing fluctuating spins 
and then the FR angle is measured without echo control (the free- 
induction decay case), (b) Modulation function in the free-induction 
decay case, (c) The Faraday rotation echo scheme. A linearly polarized laser 
pulse passes through a sample containing fluctuating spins. After 
interacting with the sample, the laser pulse passes through a linear 
birefringence crystal and then is reflected. The birefringence crystal flips 
the Faraday rotation angle of the laser. Then the laser pulse interacts with 
the sample again before its polarization is measured, (d) Modulation 
function /echo in the case of echo control. 



Faraday rotation echo spectroscopy (FRES). The idea of FRES is 
explained as follows. First we let the laser pulse go through the sample 
and accumulate an FR angle (0) by interacting with the spins in the 
sample. After the laser pulse interacting with the sample, it is 
reflected and passes through a linear birefringence crystaP\ The 
birefringence crystal has different refractive indices (n^ and n^) for 
the linearly polarized light with polarization parallel and 
perpendicular to the optical axis of the crystal. The original laser 
polarization (before interacting with the sample) is set to be 
parallel or perpendicular to the optical axis of the crystal. If the 
thickness of the crystal (L^^) is chosen such that the two 
orthogonal polarizations accumulate a phase difference 71, namely. 



2n(no 



-ne)—^=K, 



(3) 



where X is the wavelength of the laser in the vacuum, the polarization 
of the laser would be flipped after propagating through the 
birefringence crystal and the Faraday rotation accumulated before 
would change its sign (0^ — 0) [see Figure 1 (c-d)]. This is a 
straightforward analogue to the spin flip in spin echo. Then, we let 
the laser pulse go through the sample again and interact with the 
spins once more [see Figure 1 (c)] . As compared to the spin echo, this 
Faraday rotation echo can have an additional time delay t between 
the interactions with the sample before and after interaction with the 
birefringence crystal. By defining the modulation function [Figure 1 
(d)] 
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-1 ioYte[0,tp], 
fecho(t,T)={ +1 iorte[tp + T,2tp + T], 
0 else, 

the FR angle of the final output laser reads 



(4) 



Oechoi'^tp) 



1 

= (^t 'Ls 

Jo 



M,it)fecUt.T)dt. 



If the magnetization fluctuation of the spin system is static, the net FR 
angle after the echo would be zero. Therefore, the final FR measures 
the dynamical fluctuation of the magnetization, which is related to 
quantum transitions in the spin system. 

We define/FiD(^,T) = 1 for tE[0,tp] [see Figure 1 (b)] as the modu- 
lation function for the case of no echo control (similar to the 
free-induction decay in magnetic resonance spectroscopy). The 
expressions of the FR angle in the two different cases (echo and 
free-induction decay) are unified as 

e,(T,) = at-'Ls M,(t)f,(t,T)dt, (6) 
Jo 

where ?] = "FID", or "echo" indicates the free-induction decay or 
echo configurations, respectively. is the total time that the laser 
interacts with the spin system (i.e., TpiD = tp and Tgcho ^ 2/^). The 
fluctuation of the FR angle is determined by the magnetic noise 
spectrum as^^ 

{^el{T^))=OLh-^Ll [ [ dhdt2C,(h-t2)MtuT)frj(t2,T) 



1^ 



Q(0)Mtu^)Mt2,T)dtidt2 + a^ 



'til 2^ 



Sq{co} 



(7) 



where the filter function co ^\F^ (co) | ^ is the Fourier transform of the 
modulation function, with 



(8) 



The expression of (AOl^^^) in equation (7) is similar to the phase 
fluctuation of spins in magnetic resonance spectroscopy. The first 
term in equation (7) comes from the thermal fluctuation [Sth(co)], 
which vanishes in the case of FR echo since Q(0) is a constant. 
Therefore, the FR echo measures the quantum fluctuation of the 
magnetization. The tunable delay time t adds extra flexibility for 
studying the quantum fluctuations, as compared with the conven- 
tional spin echo. For example, a large delay time can be used to single 
out the effect of low energy excitations. 

The fluctuation of the FR angle will result in depolarization of the 
laser pulse. The degree of polarization of the laser pulse after the 
interaction with the sample is 



In 



-In 



Imax Im 



(9) 



where I^y^ax ^ (cos^A^^) and J^/„ oc (sin^A^^) are the maximum and 
minimum intensities of the light passes through a linear polarizer, 
respectively. Therefore, the depolarization is directly related to the 
fluctuations of the FR angle. 

FRES of LiHoF4. The lithium rare-earth tetrafluorides are a family of 
compounds used as a testing ground for the physics of spin models. 
All of these compounds are optically transparent, which makes them 
ideal for optical studies^". We choose the lithium holmium 



tetrafluoride crystal LiHoF4 



' as our model system in this study. 



LiHoF4 has a scheelite structure [Figure 2 (a)] with the lattice 
constants^" a = 5.175 A and c= 10.75 A. The magnetic properties 



of this compound come from the Holmium ions (Ho^^), which 
can be treated as a system of spin- 8. The interaction of the Ho^^ ions 
with the surrounding Li^ and F' ions are described by a crystal-field 
Hamiltonian Hcf (see Methods). The crystal field produces an 
energy level splitting of the Holmium ions. The Holmium spin has 
a ground state doublet and 15 excited states. The lowest excited 
energy level lies about 1 1 Kelvin above the ground state. When the 
temperature is much lower than 11 Kelvin, only the ground state 
doublets are non-negligibly populated. According to the expression 
of the Steven operators (see Methods), the dominant term in the 
crystal field Hamiltonian is ^ — (/f ) (the 2- direction is along the 
c axis of the crystal), which means that the ground state doublet are 
basically |8, — 8) and |8, + 8). The flipping between these two states 
needs to go through the excited states. Therefore, at temperature 
much lower than 1 1 Kelvin, the spins effectively have the Ising type 
of interaction. Quantum fluctuations of the spins can be induced by 
applying a transverse magnetic field along the x-axis. 
The full Hamiltonian of LiHoF. is"''"''"' 



H= J2 [HcFOi) + Ah-li-g^Bh-^] 



(10) 



where J/ and 1/ are the electron spin (/= 8 with^= 5/4) and the nuclear 
spin {1=712) of the /-th ion, respectively. The magnetic interactions 
in LiHoF4 include the long-range dipole interaction between the 

Holmium magnetic moments, with [Dij]^^ = D^piij) being the dipole 
sum, the exchange interaction (J12 = 0.6 fieV ~ 0.007 Kelvin) 
between Holmium spins in the nearest neighbors, and the isotropic 
hyperfine interaction between the nuclear and electron magnetic 
moments on the same site {A = 3.36 /neY ~ 0.04 Kelvin). In the 
case of zero field, LiHoF4 forms a ferromagnetic order at the critical 
temperature 1.53 Kelvin. When a transverse field is applied, the 
ordering temperature gradually approaches to zero at the critical 
field Bc = 5 Tesla. In Figure 2 (b), the expectation value of the single 
site angular momentum operator /f is calculated with the mean-field 
approximation using the full Hamiltonian in equation (10). The 
appearance of the spontaneous magnetization signatures the emer- 
gence of the ferromagnetic order. The paramagnetic-ferromagnetic 
phase boundary agrees with the phase diagram obtained in previous 
works (see, e.g., Ref. 31). Figure 2 (c-e) shows the mean field results 
of the correlation function of the magnetic fluctuations (A/^(t) A/f ) 
as a function of the transverse field B^ and the time t for various 
temperatures. In the long time limit (in which the time is longer than 
the inverse energy gap, here ^ > 0. 1 ns in our estimation), the onset of 
oscillations in the paramagnetic phase indicates the phase boundary 
in Figure 2 (b). Note that in principle, the long time limit means that 
the interacting time needs to be longer than the inverse-interaction 
energy in the spin system for the critical features to be revealed. But in 
the LiHoF4 system, the hyperfine coupling between the electron spins 
and nuclear spins causes an energy gap £gap at the phase boundary. In 
this case, the critical features appear when t > H/E^^^. It should be 
noted that in experiments^^, a finite gap is observed but with a dif- 
ferent value compared to the one in this paper. Calculation beyond 
the mean-field approximation is needed to determine more accur- 
ately the energy gap. 

The FR echo signal is evaluated according to equation (7). The 
correlation function Q is obtained with the mean-field approxi- 
mation using the full Hamiltonian, which is related to the single spin 
correlation functions by 

Q{t) = (AM,(0)AM,(0> =N-'{Ajf{0)Ajf{t)). (11) 

Note that in the mean-field approximation, the correlations between 
spins at different sites vanish (i.e., (A/f (0)A/^^(t)\ = 0 for / 7^ j). 
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5,(Tesla) 

Figure 2 | Phase diagram of LiHoF4 and magnetic correlation functions, (a) Lattice structure of LiHoF4. (b) Magnetization of LiHoF4 as a function of 
temperature and transverse magnetic field. Sudden appearance of magnetization along the z-axis indicates the formation of the ferromagnetic order. 
The correlation function {Ar(t) AD is plotted as a function of the transverse field and time tat temperature (c) 0.3 Kelvin, (d) 0.9 Kelvin and (d) 1.35 
Kelvin. The corresponding critical fields are indicated as open circles in (b) and dashed lines in (c) - (e). An obvious oscillation feature appears in the 
paramagnetic phase. 



Figure 3 (a) plots the FR angle fluctuation under the free-induction 
decay configuration as a function of the temperature and the trans- 
verse field. When the system is in the paramagnetic phase, the fluc- 
tuation of the magnetization is suppressed by the transverse field and 
is less sensitive to B^. A critical feature, namely, the sudden change of 
slope [see inset of Figure 3 (a)], appears at the phase boundary. 
Figure 3 (b) shows the FR echo signal (for delay time t = 0). After 
the thermal fluctuation is removed by the echo control, a peak feature 
is observed at the phase boundary, where the quantum fluctuation 



diverges [see inset of Fig 3 (b)]. According to the estimation (see 
Methods), for a laser with wavelength X = 435 nm and cross-section 
area A = 10 |im^ the strength of the echo signal is about 10"^ 
(degree)^, which is experimentally observable (see e.g., Ref 32). It 
should be pointed out that the mean-field approximation used in the 
calculation underestimates the fluctuation near the phase boundary 
since the long-range correlations of the fluctuations emerge at the 
phase transitions but the mean- field approximation considers only 
the local correlations. 




Figure 3 | Fluctuations of the Faraday rotation angle, (a) The FR angle fluctuation (A^pj^) in the free-induction decay configuration (no echo control) is 
plotted as a function of temperature and transverse field, for interaction time rpiD= 0.3 ns. The inset shows (A^pj^) as a function of the transverse field at 
temperature 0.9 Kelvin, (b) The fluctuation of the FR angle under echo control (A^^^j^^) as a function of temperature and transverse field, for 
interaction time Techo—O-^ ns and delay between two interaction intervals t = 0. The inset shows (A^^^^j^^) as a function of the transverse field at 
temperature 0.9 Kelvin. The sample thickness Ls = 0.5 cm, and laser cross sectional area A = 10 |^m^ are assumed in all the estimations. The white curves 
in (a) and (b) indicate the phase boundary. 
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Figure 4 | Noise spectra and the effect of different interacting times, (a) Noise spectra Sq at temperature 0.9 Kelvin for various transverse magnetic fields, 
(b) The lowest excited state energy as a function of the transverse field at temperature 0.9 Kelvin. The purple circles in (a) and (b) mark the lowest energy 
excitation at transverse field around = 3.2 Tesla. (c) (A^^^j^^) as a fianction of transverse field at temperature 0.9 Kelvin for various interaction time 
Techo and delay time t = 0 ns. The sample thickness is set as L5 = 0.5 cm, and laser cross sectional area is set as A = 10 |am^. 



Discussion 

To better understand the relation between the phase transitions and 
the features in the FRES, we study the magnetic noise spectra of the 
LiHoF4 system. Figure 4 (a) shows some typical examples of the noise 
spectra at temperature 0.9 Kelvin. The gap of the system reaches the 
minimal value (—6.7 ^eV) around = 3.2 Tesla [see Figure 4 (b)], 
which indicates the onset of the phase transition. Therefore, the low 
frequency components of Sq have larger contribution to the Faraday 
rotation fluctuation when the external field is closer to the critical 
point. Since we are considering the case of finite temperature, there 
also exist some lower frequency components of noise spectra due to 
the transitions between the excited states. But the contributions from 
those components are very small due to the small distribution prob- 
abilities of the excited states. Therefore, the critical features around 
the phase transition point are mainly determined by the transition 
between the ground state and the lowest excited state [as indicated by 
the purple dashed circle in Figure 4 (a-b)] , when the interacting time 
is comparable to or longer than the inverse excitation gap (the "long- 
time" limit). 

Figure 4 (c) shows the dependence of the peak features of (AOl^^^) 
on the total interaction time Techo between the laser pulse and sam- 
ple. When Techo is long enough, a peak feature emerges at the phase 
transition point. The peak becomes sharper as the interacting time 
^echo further increases. The interaction time T^ = 0.3 ns in Figure 3 
(a-b) is chosen long enough for the critical features to pronounce. 

As compared with the conventional spin echo in magnetic res- 
onance spectroscopy, the FRES has extra controllability to engineer 
the filter function by tuning the delay time between the interaction 
intervals. To show this controllability, we study the effect of the delay 
time. Figure 5 (a) plots the FR echo signal with a delay of 0.25 ns. The 
effect of the delay time could be considered as expanding the modu- 
lation function to a longer time interval and making the low fre- 
quency excitations more important [see Figure 5 (b)]. Therefore, 
the critical feature at the phase boundary persists as we increase 
the delay time [see Figure 5 (c)]. In addition, a long delay time can 
also reveal the critical features around the phase transition points. In 



Figure 5 (d), the FR echo signal of Techo — 0.1 ns is plotted with 
different delay times. By prolonging the delay time, a sharp peak 
emerges at the phase transition point. 

The FRES scheme can be straightforwardly extended to more 
complicated configurations, corresponding to different kinds of 
dynamical decoupling sequences^^'^^, by letting the laser pulse inter- 
act with the sample and the birefringence crystal multiple times. 

In this paper we have assumed a square shape for the laser pulse. In 
general, by shaping the laser pulse one can realize more complicated 
modulation functions (which would be particularly useful for, e.g., 
spectroscopy of the spin noises in the bath^^'^^). 

Methods 

Faraday Rotation. For simplicity, we consider the FR of a linearly polarized laser 
weakly coupled to a single spin J. The Faraday rotation is originated from the relative 
phase shift of the two circular polarized mode (with photon annihilation operators b+ 
and b-). The effect Hamiltonian is 

v=xnbib+-bib_), (12) 

where x <^ i^tp)~^ is the coupling strength. The weak coupling condition is X^p-^^- 
The initial state of the whole system is on a factorized state 

|<^o>=e^'«-^*'^|0>(x)|/,m>, (13) 

where bH = {b+ +b-)/ a/2 is the annihilation operator of the linear horizontal mode, 
that is, the laser is initially on a linearly polarized coherent state, and \J,m) is an 
eigenstate of with eigenvalue m. The evolution of the state is 

\(j>{t))=e''^'^mm-^-'-\Q)®\J,m), (14) 

where be^ = cos dmbn — sin Ombv is the annihilation operator of the linearly 
polarized mode with the polarization rotated by an angle 9m{t) = x^^- Here, 

by = -j^ (b+ —b-)is the annihilation operator of the vertically polarized mode. The 

ensemble measurement of the FR angle gives 

{0(t))=x{lz)t. (15) 
When a linearly polarized laser couples to N fluctuating spins with Hamiltonian 
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as a function of temperature and transverse field, for interaction time Techo^O-^ ns and delay 
0.25 ns. The white curve indicates the phase boundary, (b) The filter function w 



^l^echol^ for '?'echo=0-3 ns and 



various delay times t. (c) (A^^^^j^^) as a function of transverse field at temperature 0.9 Kelvin for Techo =0-3 ns and various delay times t. (d) (A^^^j 



function of transverse field at temperature 0.9 Kelvin for T^cho =0.1 ns and various delay times t. The sample thickness is set as L5 = 
and Ls = 0.1 cm in (d). The laser cross sectional area is set as A = 10 |am^ in all the figures. 



0.5 cmin(a) &(c), 



V{t) = xNMMib^^b+ -blb_), (16) 
the Faraday rotation angle is 

eit) = xNt M,(t')dt', (17) 
Jo 

where Mz(t) = N~^ ^ Jf(t) is the magnetization. Considering x oc {Atp)~^ and AT oc 
ALs, we obtain equation (1). 

Crystal-field Hamiltonian. The crystal-field Hamiltonian is^^ 

Hcp= 5?0?+ ^5/(c)0?(c) + 5^(5)0^(5), (18) 

/ = 2,4,6 1 = 4,6 

with the Steven's operators OJ" and the corresponding crystal-field parameters B'p. 
The S4 symmetry of the Ho -ion surroundings is responsible for the terms in equation 
(18). The specific forms of the Steven's operators are listed below. 

02 = 3(/f)'-/(/+l), 

O2 = 35(/f)'-30/(/+l)(/f)'+25(/f)'-6/(/+l) + 3/2(/ + l)^ 
0« = 23l(/f)'-315/(/+l)(/f)' + 735(/f)' + 105/2(7+ l)2(/f)' 

-525/(7+ l)(7f)'+294(/f)' -5/^(7+1)' +40/2(7+1)' -60/(7+1), 

oi(c)=i[(/^.-)V(/r)^], 

Ot(c)=l[{jrT + {jrY] [ll(/f)^-/(/+l)-38]+H.c., 
Ot(s)=^.[{jn'-{hn'] [ll(/f)^-/(/+l)-38]+H.c., 
where Jr=Jf±iJl 

In this paper, we use the same crystal-field parameters as in Ref. 28, as listed below. 
bI=-60 fieV, bI = 0.35 fieV, 5^ = 3.6 fieV, 
bI = 0.0004 fieV, Blic) = 0.07 fieV, 5^(s) = + 0.0098 fieV. 



Mean field calculation. The results we present in Figs. 2-5 are calculated with the 
mean field theory. The full Hamiltonian of the LiHoF4 system in equation (10) 
includes two parts. The non-interaction part 



which is a single-ion Hamiltonian, and the interaction part 



(19) 



(20) 



which is treated with standard mean-field approach (see, e.g., Ref 35). The mean-field 
Hamiltonian is obtained by neglecting the two-site fluctuation terms 



V 

-^E^>2[2(J;)-Ji-a)-(J/)i- 



(21) 



After replacing Hi with H^^, the Hamiltonian of the LiHoF4 system becomes a 
single-ion Hamiltonian, which is then exactly solved by numerical diagonalization. 
The mean field (J^) is obtained by iteration. 



Estimation of the signal strength. The coupling strength a is evaluated with 



{M,)Ls' 



(22) 



for static magnetization. According to Ref 20, 9/Ls ~ 9260 degree/cm in the case of 
saturated polarization ((M^) = 8) for a laser with wavelength 435 nm. Therefore a = 
1157.5 degree/cm is used for signal strength estimation in this paper. 

Since the correlations between different electron spins are neglected in the mean- 
field approximation, we have Cz(0 = N~^(A/f (0)/f (f)). Therefore, 

= ^ [ ('^^ {Aj\h)Aj\t2))fn{tuT)m2,T)dhdt2, (23) 
P^'^p J Jo 

where p ~ 1.39 X 10^*^ m"^ is the density of spins in LiHoF4 system. The results in 
Figures 3-5 are estimated according to the above equation. It should be noted that, the 
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mean-field approximation is only valid away from the critical points. In this paper, 
the mean-field results are used to illustrate the qualitative properties of the FR echo 
signals around the phase boundary. 
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